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Introduction 

Diffusion-limited aggregation (DLA) is a model for dendritic growth in 
which a particle executes a random walk in the lattice Z d beginning "at 
infinity" and ending when first adjacent to a region AcZf The terminus 
of the random walk is then adjoined to the region A, and the procedure is 
iterated. It is conjectured |2| that the resulting region A has fractal dimension 
exceeding 1. Internal DLA (IDLA) is a variant in which the random walks 
begin at some fixed point in A and end when they first leave the A. Unlike 
DLA, internal DLA does not give rise to fractal growth. Lawler et al. (|17j. 
1992) showed that after n walks have been executed, the region A, rescaled 
by a factor of n l ' d , approaches a Euclidean ball in R d as n — > oo. Lawler 
([T5j. 1995) estimated the rate of convergence. 

We propose the following deterministic analogue of IDLA. First, a cyclic 
ordering is specified of the 2d cardinal directions in Z d . A "rotor" pointing in 
one of the 2d cardinal directions is associated to each point in A. A particle 
is placed at the origin and is successively routed in the direction of the rotor 
at each point it visits until it leaves the region A. Moreover, every time the 
particle is routed away from a given point, the direction of the rotor at that 
point is incremented by one step in the cyclic ordering. When the particle 
reaches a point not in A, that point is adjoined to the region and given the 
first rotor direction in the ordering. 

This "rotor-router" growth model is similar in many ways to Engel's 
"probabilistic abacus" ([IB], 1975) and to the abelian sandpile, or "chip- 
firing," model introduced by Bak et al. ([T], 1988) and studied by Dhar ([§], 
1990) and Bjorner, et al. ( 5J, 1991). In the abelian sandpile model, the 
points of A are labeled by nonnegative integers, considered as representing 
the number of grains of sand at each point. A single grain of sand is placed 
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at the origin, and at each step, every point occupied by at least 2d grains of 
sand ejects a single grain to each of its 2d neighboring lattice points. Once 
the system has equilibrated, with every site occupied by at most 2d — l grains 
of sand, the process is repeated. Sandpiles have been studied in part from the 
point of view of complex systems, where the interest lies in the "self-organized 
criticality" of the model (cf. PQEHH]). In contrast to diffusion-limited aggre- 
gation, however, many of the fundamental mathematical properties of which 
remain conjectural, sandpiles are tractable. 

Much of the mathematics concerning the sandpile model rests on the 
abelian nature of the model. This refers to the nontrivial fact that if grains 
of sand are deposited in turn at two different points, allowing the system to 
equilibrate both before and after the second grain is deposited, the resulting 
configuration does not depend on the order in which the two grains were de- 
posited. It is a consequence of a general result of Diaconis and Fulton 
1991) that the rotor-router model has the same abelian property. Another 
result of PI] implies that internal DLA is abelian in the sense that the prob- 
ability that two random walks with different starting points will terminate 
at a pair of points x, y is independent of the order in which the walks are 
performed. 

Despite this common abelian property, there are substantive differences 
between IDLA and the sandpile model. While no asymptotic results are 
known for the sandpile model, it appears likely that its asymptotics are 
not spherical. In two dimensions, for example, numerical data indicate that 
sandpiles may have polygonal asymptotics. This is to be contrasted with the 
main result of [17|, in which IDLA is shown to have spherical asymptotics. 
The rotor-router model may bridge the gap between sandpiles and IDLA. Like 
sandpiles, the rotor-router model is deterministic; but we conjecture that the 
asymptotics of the rotor-router model, like those of IDLA, are spherical. 

This paper is intended to serve, first, as a thorough introduction to the 
rotor-router model, and second, as a source of conjectures and open problems 
which, it is hoped, will inspire future work on the model. Analogies with 
IDLA and with the sandpile model are emphasized throughout. The paper 
is structured as follows. 

In section 1, we outline some preliminary definitions and prove an impor- 
tant finiteness lemma. 

In section 2, we study the one-dimensional (d — 1) rotor- router model in 
some additional generality. Fixing positive integers r and s, when a particle 
leaves the interval A through its left endpoint, A is extended by not just 
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one but r sites to its left; and when a particle leaves A through its right 
endpoint, A is extended by s sites to the right. We reduce to a set of states 
describing the eventual behavior of the system. By identifying this collection 
of states with a suitable subset of the integer lattice Z 3 , we show that the 
process of adding a particle at the origin and allowing it to equilibrate can 
be realized by a piecewise linear function on Z 3 . We find a set of invariants 
of this function that is sufficient to distinguish between all its orbits. Our 
main invariant closely resembles an invariant of the continuum limit of IDLA 
in one dimension. 

In section 3, we analyze the limiting behavior of the interval A in the 
generalized one-dimensional model. Propp (2001) conjectured that after t 
iterations, if A is the interval [x(t),y(t)], then the quantity /i(t) = x(t)y/s + 
y{t)\/r is bounded independent of t. We prove this, and show furthermore 
that the limit points of the sequence \i are confined to an interval of length 
r^fs + Sy/r, which is best-possible given that x and y change by increments 
of r and s, respectively. We show that in a special case, the n-th particle 
added ends up on the left or right of the interval accordingly as the n-th 
term of a certain Sturmian sequence with quadratic irrational slope is or 
1. Finally, we give a connection with Pythagorean triples. Given a positive 
Pythagorean triple a 2 + b 2 = c 2 , Propp observed that in the continuum limit 
of IDLA, if initially A is the interval [—a, 0], then A is the interval [— c, b] at 
time a+b—c, and conjectured that the same would be true of the rotor-router 
model in the case r = s = 1. We prove this conjecture. 

In section 4, we consider the higher-dimensional lattices Z d , d > 2. We 
give an estimate for the center of mass of the region A and prove a weak 
form of the conjecture that the limiting shape of the rotor-router model in 
two dimensions is a disc. 

In section 5, we discuss some outstanding conjectures pertaining to the 
one-, two- and three-dimensional models. 

1 Preliminaries 

We denote by Z the set of integers, N the nonnegative integers, and 1R the real 
numbers. If x < y G Z, we denote by [x, y] the interval {z G Z : x < z < y}; 
we adopt the convention that [x, y] is empty when x > y. 

Let ei, e d be the standard basis vectors for the integer lattice Z d , and 
denote by E d = {±e 1: . . . , ±e d } the set of cardinal directions in Z d . Suppose 
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we are given an arbitrary total ordering < of the set E^. write Ed = {e^}^ 1 
with e < • • • < e 2d -i- 

Let A C be a finite connected region of lattice points containing the 
origin. A state of the rotor-router automaton can be described by a pair (x, I), 
where x G A represents the position of the particle, and I : A — > [0, 2d — 1] 
indicates the direction of the rotor at each point. Define g(x, I) = {x+eu x \, l x ), 
where l x is the labeling 



l(x') + 1 (mod 2n) if x' = x; 
l(x'), else. 



This is the state given by routing the particle in the direction of the rotor e^ x ) 
and then changing the direction of the rotor. Write g n (0,l) = (x n ,l n ). The 
sequence of points (x , Xi, . . . ) is a lattice path in A, possibly self- intersecting, 
beginning at the origin. We denote this path by p = p(A,l). Lemma ll.ll 
below, shows that the path p leaves the region A in finitely many steps. Let 
iV be minimal such that xn ^ A. We define f(A, I) = (A U {xn}, I')-, where 



l\x) 



In(x) if x E A; 
0, if x = xn- 



Then / describes the entire process of adding a particle at the origin and 
allowing the system to equilibrate. 

Lemma 1.1. The lattice path p(A,l) leaves the region A in finitely many 
steps. 

Proof. If not, the path would visit some point x G A infinitely many times; 
but then it would be routed infinitely many times to each neighbor of x. 
Inducting along a path from x to a point outside A, we conclude the path 
does after all leave the region after finitely many steps. □ 



2 Rotor-router dynamics in one dimension 

In one dimension, rotors alternate between the two directions left and right; 
we denote these by L and R, respectively. We introduce the following gen- 
eralization of the rotor-router automaton in one dimension. Let r and s be 
positive integers. When the lattice path reaches an unoccupied site x < 0, 
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all r sites in the interval [x — r + 1, x] become occupied; similarly, if the path 
reaches an unoccupied site y > 0, the s sites [y, y + s — 1] become occupied. 
In either case, the newly occupied sites are initially labeled R. 

Given integers x < and y > 0, we denote by y) the set of all states 
of the automaton for which the set of occupied sites is the interval [x,y]. 
Then y) is naturally identified with the set of maps [x,y] — > {i?, L}. 
Let X = |J a . <0<?/ y), and denote by / = / ns : S — > X the map on states 
given by adding a single particle at the origin and allowing the system to 
equilibrate. 

2.1 Recurrent states 

Lemma 2.1. Let a G y) fre any state. There exists N e N snc/j £/ierf 
/^>)eE(^yOW<z ; y'>y. 

Proof. Let M = 1 + max(|x|,y), and let TV = 2 M . After iV particles have 
been deposited and allowed to equilibrate in turn, the origin will have been 
visited a total of at least N times. Since the rotors at each site alternate 
pointing left and right, it follows by induction on \k\ that if < M, then 
the site k will be visited at least 2 M ~ k times. Taking k = ±M proves the 
lemma. □ 

Given a state a e y), let > u± > ■ ■ ■ > u m be the sites to the left 
of the origin labeled R, and let < V\ < ■ • • < v n be the sites to the right 
of the origin labeled L. Additionally, define u = v = 0, u m+ i = x — 1, and 
v n+ i = y + 1. Then the path p{a) can be described as follows. 

Lemma 2.2. // cr(0) = L, then the path p(a) travels directly left from the 
origin to U\, then right to v±, left to u 2 , right to v 2 , and so on, until it reaches 
either u m+1 or v n+1 , at which point it stops. If <j(0) = R, the same is true 
interchanging left with right and Ui with v^. In particular, if m < n, the path 
will terminate at u m+ \; if m > n, it will terminate at v n+ ±; and ifm — n, it 
will terminate at u m +i or v n+ i accordingly as cr(0) = L or cr(0) = R. 

Proof. Suppose er(0) = L. Induct on k to show that when the path first 
reaches the site Vk, it travels left from f fc to then right from to 

i>fc+i. When the path first arrives at Vk, it must have previously reached 
Vk-i, so by the inductive hypothesis, the path has come directly right from 
u k to Vk, hence the entire interval [iifc,i>fc] is now labeled L. Since also, by 
definition, the interval [u k+1 + l,u k — 1] is labeled L, the path now travels 
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directly to the left until it reaches the site u^+i- Now the interval [uk+i, Vk] is 
entirely labeled R, as is the interval [vk + 1, v k+1 — 1] by definition; so the path 
travels directly right from Uk+i to v^+i, and the inductive step is complete. 

The proof in the case <r(0) = R is identical, interchanging the roles of Uk 
and Vk- □ 

Remark. It follows that the N in Lemma f2. II can be taken substantially less 
than 2 M . If, for example, the path p{f l {o~)) terminates to the right of the 
origin for i = 0, . . . , k, then for i > 1 we have n (f l (a)) = y + (i — l)s, 
and so m(/ I+1 (cr)) = m(f l (o~)) — y — (i — l)s. Thus, to ensure that at 
least one particle terminates on the left it suffices to take N large enough 
so that Ny + s jV ( j ^~ 1 ) > \x\. Likewise, to ensure that at least one particle 
terminates on the right it suffices to have N\x\ + r Af ( J ^ +1 ) > y Certainly 
N = 1 + A/2M/min(r, s) is enough. 

Let Rec(x, y) C S(x, y + s — 1) be the set of states for which there exists 
an integer < i < y — x + 1 such that the interval [x, x + i — 1] is entirely 
labeled R, the interval [x + i, y — 1] is entirely labeled L, and the interval 
[y,y+s — 1] is entirely labeled R. Putting j = y — x — i, we say as a shorthand 
that these states are of the form R l DR s for nonnegative integers i and j. 
Let Rec = \J x < < y Rec(x,y). 

We are interested primarily in the eventual behavior of f T ^ s as a dynamical 
system on S. The following proposition shows that for these purposes, it 
is sufficient to consider the states in Rec. These states will be called the 
recurrent states of f r>s . 

Proposition 2.3. The set Rec is closed under f r ^ s . Moreover, for any state 
cr £ S, there is an integer N > such that f^ s (cr) G Rec. 

Proof. Let a G Rec(x, y). Suppose first that the path p(a) terminates at 
y + 1. Then by Lemma 12.21 p(o~) will reach a site u < and travel directly 
right from there to the unoccupied site y + 1; in particular, the path never 
visits the interval [x,u — 1], and f(a) retains its original labels from a on 
this interval. Since u is in the R 1 block of a, all of these labels are R. Also, 
because p(a) travels to the right from u to y + 1, the interval [it, y] is entirely 
labeled L; and the s newly occupied sites in the interval [y, y + s — 1] are 
labeled R. Thus / r , s (cr) has the form R l DR s , as desired. On the other hand, 
if p(cr) terminates at x — 1, then by Lemma 12.21 it will reach some v > 
and travel directly left to the unoccupied site x — 1. In this case, the interval 
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[x — r, v] is entirely labeled R, while the sites in the interval [v n + 1, y] retain 
their original labels from a. Since a G Rec, the labels on this interval are of 
the form L^R S , so again f r . s {o~) has the desired form. 

It remains to show that for any o G X, some iterate of / takes a into Rec. 
By Lemma 12.11 there exists N such that p (f N ~ 2 ((j)) terminates on the left 
and p (/ Ar_1 (c)) terminates on the right. Since p (f N ~ 2 {o~)) terminates on 
the left, we have f N ~ 1 (a)(t) = R for t < 0. Since p(/ 7V_1 (cr)) terminates 
on the right, there is some point U{ < such that f N (u) retains its original 
labels from f N ~ 1 (cr) on the interval [x, Ui — 1] and has the form D R s on the 
interval [ui,y]; hence f N {c) is of the form R l UR s as desired. □ 

2.2 A piecewise linear function on Z 3 

A state a G Rec has the form R l L 3 R s , and so is determined by the pair of 
integers i and j (recall that s was fixed at the outset, independent of a). If we 
intend to compute / r>s (o"), however, then the origin must be distinguished so 
that we know where to initiate the path p(cr). With the origin distinguished, 
Rec becomes a three-parameter family of states. There are several reason- 
able parameterizations of Rec, but the one that simplifies computation most 
effectively is to let x, y, z G Z be the first occupied site, the first site in the 
final R s block, and the first site in the V block, respectively (if j = 0, we 
adopt the convention that z = y). In this way, we identify Rec with the set 
of integer triples (x, y, z) G Z 3 satisfying x < 0, y > and x < z < y. Our 
next proposition determines / rjS explicitly as a piecewise linear function on 
these triples. 



Proposition 2.4. f T)S is given on Rec by 
f rjS (x,y,z) = 



y + s, z — y) if x + y < z. 
x — r, y, z — x + 1) if x + y > z, 



Proof. Consider first the case z > 0, i.e. the case when the origin is initially 
labeled R. Then with m and n defined as in Lemma 12.21 we have m = —x 
and n = y — z. By Lemma \'2. 21 the path p = p(x, y, z) ends up on the right if 
and only if m > n, or x + y < z. In this case, once the path reaches the site 
u n = —n = z — y, it travels right until it reaches an unoccupied site. The 
first site labeled L is then z — y, so 

f(x,y,z) = (x, y + s, z- y). 
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On the other hand, if x + y > z, then the path travels left directly from 
the site t> m +i = z + m until it reaches an unoccupied site. If f m +i < V~ 1, the 
first site labeled L is then v m +i + 1 — z + m + 1 — z — x+1; if v m+ i = y — 1, 
then there are no sites labeled L, so our convention dictates that z' — y — 
Vm+i + 1 — z — x + 1] hence 

f(x, V,z) = (x-r, y, z-x + 1). 

It remains to consider the case z < 0. In this case, m = z — x and 
n = y — 1. The path ends up on the right if and only if m > n, or x + y < z, 
as before. If it ends up on the right, the path travels directly right from 
u n+ i = z — n — 1 = z — y to the unoccupied site, so the first site labeled L is 
z — y, as desired. If the path ends up on the left, then it travels directly left 
from v m = m = z — x to an unoccupied site. In the case that v m < y — 1, the 
first site labeled L is then z — x + if v m = y — 1, there are no sites labeled 
L, and by convention z' = y = z — x + 1. This completes the proof. □ 

We will denote the states (x,y + s, z — y) and (x — r, y, z — x + 1) by 
f + (x,y,z) and f~(x,y,z), respectively. Thus f(x,y,z) = f + (x,y,z) or 
f~ (x, y, z) accordingly a.s x + y<z or x + y>z. 



2.3 An invariant 

Consider an analogous generalization of stochastic IDLA in one dimension, 
in which r or s sites become occupied accordingly as the random walks ter- 
minate on the left or right sight of the interval. Each random walk is a 
"gambler's ruin" problem (see, for example, [13]), terminating on the right 
with probability tM^, where [x,y] is the interval of occupied sites. Thus the 

limiting value of the ratio ^ as time goes to infinity satisfies 

\x\ 



ry/(\x\ 


+ y) 




x\/(\x 


\ + y) 



hence 

\x\ ir 



as t — ► oo. 

y v s 

In fact, in the continuum limit of IDLA, as the frequency with which particles 
are dropped is taken to infinity and the interval is rescaled appropriately, the 
model becomes deterministic and the quantity sx 2 — ry 2 is exactly conserved. 



8 



This suggests that the quantity sx 2 — ry 2 is likely to be close to invariant in 
the rotor-router model as well. As the following lemma shows, this is indeed 
the case. 

Lemma 2.5. The function 

g r , s (x, y, z) = sx 2 — ry 2 + (r — 2)sx + rsy — 2rsz (1) 

is invariant under f rjS . 

Proof. Compute 

9(f~(x,y,z))-g(x,y,z) = g(x - r, y, z - x + 1) - g(x, y, z) 

= s(x — r) 2 — sx 2 — (r — 2)sr + 2rs(x — 1) 
= — 2rsx + r 2 s — r 2 s + 2rs + 2rsx — 2rs 
= 0; 

g(f + (x,y, z)) -g(x,y, z) = g(x, y + s, z - y) - g(x,y, z) 

= 0. □ 



The following section is devoted to showing (Theorem 12.8(1 that this in- 
variant, together with the congruence classes of x (mod r) and y (mod s), is 
sufficient to distinguish between all orbits of /. Later, in sections 3.1-3.3, we 
use this invariant to derive a variety of bounds on the growth of the interval 

2.4 Classification of Orbits 

Since g is linear in z, it follows from Lemma (2.51 that given any x,y,n G Z, 
there is at most one state a = (x,y,z) G Rec for which g(a) = n; namely, 
set 

sx 2 — ry 2 + (r — 2)sx + rsy — n 
z = z n (x,y) = — 

x(x + r — 2) y(y — s) n 
2r 2s 2rs 

if this is an integer and satisfies x < z < y. 
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Lemma 2.6. Let F(x,y) = z n (x,y) — x and G(x,y) = z n (x,y) — y. Then 
for all x < 0, y > we have 

(*) F(x,y) < F(x - r,y); 
(ii) G(x,y) < G(x - r,y). 
(Hi) F(x,y) > F(x,y + s); 

Proof. From (J2J) we have 

x(x-r-2) y(y-s) n 

r (x, y) = . 

v ,y; 2r 2s 2rs 

Thus, for fixed x, F is a quadratic in ?/ with maximum at y = |; and for 
fixed y, F is a quadratic in x with minimum at x = 1 + |. This proves (i) 
and (ii). Likewise, 

G( s = x ( x + r ~ 2 ) _ + g ) _ 

2r 2s 2rs' 

so for fixed y, G is a quadratic in x with minimum at x = 1 — |, and this 
proves (iii). □ 

We adopt the notation (x, y, n) as a shorthand for the state (x, y, z n (x, y)). 
We will say that two states (x, y, z) and (x', y', z') are congruent if x = x' 
(mod r) and y = y' (mod s). By the congruence class of (x,y,z) we will 
mean the pair of integer congruence classes (x mod r, y mod s). Trivially, 
congruence class is invariant under / rjS . Also, notice that if z n (x,y) G Z 
and (x',y',n) = (x,y,n), then z n (x',y') G Z by (0)- Given a state cr = 
(x , yo; n ) ^ R- ec an d fixing y > y congruent to y (mod s), let 

A n (y) = {x < : x = x (mod r), x < z n (x, y) < y}. 

Then (x, y, n) G Rec if and only if x G A n (y). 

If m G N and a < b G Z with a = b (mod m), we denote by [a,6] m 
the spaced interval {a, a + m, a + 2m, . . . , b}. This interval is said to have 
increment m. The following lemma describes the set of pairs (x, y) for which 
the triple (x, y, n) is in Rec. 

Lemma 2.7. v4 n (y) o nonempty spaced interval with increment r. Fur- 
thermore, writing A n (y) = [a~ (y) , a + (y)] r , the upper endpoint a + satisfies 
a + (y + s) < a + (y). 
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Proof. Notice that x G A n (y) if and only if 

F(x,y)>0>G(x,y). (2) 

By parts (i) and (ii) of lemma 12.61 the set of x < congruent to Xq (mod r) 
satisfying ((2)) is a spaced interval of increment r. This interval is nonempty 
since by lemma I2~T1 there exists iV G N such that f N (o~o) has right endpoint 

y- 

Now since a + {y) + r A n (y) and G(a + (y) +r,y) < G(a + (y),y) < 0, it 
must be that F(a + (y) +r,y) < 0. Thus, by part (hi) of lemma I2~B1 for any 
k G N we have 

F(a + (y) + kr,y + s)< F(a + (y) +r,y) < 0, 

hence a + (y + s) < a + (y) as desired. □ 

If O C Rec is an orbit of /, we write g{0) for the constant value of g 
on O, and c(O) for the common congruence class of the elements of O. Our 
next result shows that g and c are a complete set of invariants for / in the 
sense that no two orbits of / have the same pair (g, c). 

Theorem 2.8. Suppose o~, a' G Rec are congruent states. Then a and a' 
are in the same orbit of f rjS if and only if g r s ((r) = g r ^ s (a'). 

Proof. The "only if" direction is Lemma 12.51 For the "if" direction, let n = 
g r ,s(o~), and let O be the set of states (x, y,z) = o satisfying g(x,y,z) = n. 
If a = (x,y,z) G O, then certainly f(a) G O, so O contains at least one 
of the two states f~(at) = (x — r, y, z — x + 1), f + (ct) = (x,y + s, z — y) in 
the definition of the piecewise linear function of Proposition 12.41 But f + {ot) 
is a legal state if and only if z — y > x, while /~ (a) is legal if and only if 
z — x+ 1 < y, and these two conditions are mutually exclusive. Thus exactly 
one of the states f + (a) and /"(«) is in O. 
By Lemma f2. 71 for any y we can write 

A n (y) = [a~(y),a + (y)] r , 

Denote by a ± (?/) the states {a ± {y),y ) n). Since a~(y) is the lower endpoint of 
the interval A n (y), we have f~(a~(y)) ^ O, hence f + (a~(y)) = (a~(y),y + 
s, n) G O. Thus 

a-{y)<a + {y + s). (3) 
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But by lemma 12771 a + (y + s) < a + (y), so the state (3 := (a + (y + s),y,n) is 
in O. Since 

f + (f3) = ct + {y + s) G O, 

we have f~(P) £ O, and it follows that equality holds in 0. Letting cr G O 
be the state with minimal we conclude that every a G O is f k {o~o) for 

some fc. □ 



3 Bounds for the one-dimensional model 

To any state a G Rec we associate an infinite lattice path P = P rs (a) in 
the first quadrant, starting at the origin, whose n-th step is up or to the 
right accordingly as the n-th particle added to a ends on the right or left of 
the interval. This section is devoted to bounding — and, in a special case, 
determining exactly — the shape of the path P. 



3.1 Linear bounds 

The following result shows that the path P r ^ s (o~) is bounded between two 
parallel lines of slope a = In particular, for any e > 0, there is a line 
I of slope a and a translation £' = i + (1 + e, — 1 — e) of i such that all but 
finitely many steps of P lie between I and £' (figure 1). 

Theorem 3.1. Let a G Rec be any state, and set a t = (x(t),y(t), z(t)) 
— fr,s{o'{t — l)). Then l^v^ + ^v^l ^ s bounded independent oft. Specifically, 
given any e > there exists N G N such that for all t > N , 

" <x^ + yyfs<^ ^ w - . (4) 

Remark. The difference between the upper and lower bounds in (J1J) ap- 
proaches r^fs+s^ as e — > 0. This is best-possible, since the difference of the 
bounds must be at least the difference of X\fs+{y+s) ^fr and (x—r)\fs+yy/r. 

Proof. Write u — |, v = |. Let C = g rjS (a). Since g TtS is invariant, we have 
sx 2 + (r — 2)sx = ry 2 — rsy + 2rsz + C (5) 
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Figure 1: The lattice path P5 i i(0,0,0) is asymptotically bounded between two lines 
of slope y/B separated by a translation of (1, —1). It crosses the lower line at the points 
(1, 1), (5, 10), (221, 493), (1513, 3382), .... The dotted steps show the unique path bounded 
between these two lines, where it differs from P. 



at any time t. Then completing the square and using the fact that z < y, we 
obtain 

s(x + u — l) 2 — s{u — l) 2 = ry 2 — rsy + 2rsz + C 

< ry 2 + rsy + C 

= r(y + vf -rv 2 + C. (6) 

Let C = s(u — l) 2 —rv 2 + C. Lemma \'2. II ensures that x — > — oo and y — > oo 
as t —>■ oo. In particular, we can take t sufficiently large so that x < 1 — u, 
y > v , and r(y + v ) 2 > C . Then we obtain from © 

Vi(i-x-u) < VKy + ^ + c ( 7 ) 
< My + v) + , 

which gives a time-independent lower bound for X\fs + 2/\/r. 

Similarly, isolating the terms of © involving y, completing the square, 
and using the fact that z > x gives 

r(y — v) 2 — rv 2 = sx 2 + (r — 2)sx — 2rsz — C 
< sx 2 - (r + 2)sx - C 

= s(x - u-1) 2 - s(u+l) 2 -C. (8) 
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Write C" = rv 2 — s(u + l) 2 — C. By Lemma 12.11 we can take t sufficiently 
large so that x < —u — 1, y > v and s(x — u — l) 2 > C", and estimate 



y/r(y-v) < ^ s{x - u - l) 2 + C" (9) 
< yfi(u + 1 - x) + \/\C ir \, 

which gives a time-independent upper bound. 

To show (jU, choose i sufficiently large so that ^^ y+v ^ and ^ u+ [_ x ^ are 
strictly less than e. Then using first-order Taylor estimates for the square 
roots in (JJJ) and Q, we obtain 

\C'\ e 
y/s(l - x - u) < ^r(y + v) + —=- — ■ — - < y/r(y + v) + -; 

2y/r(y + v) 2 

\C"\ e 

y/r{y -v)< y/s(u + l-x) + — —. — — r < y/s(u + 1 - x) + -; 

l\J s[u + 1 — x) 2 

and hence 

2 = 0- ~ u )Vs ~ vVr - - 



< x\fs + y^/r 

l + it)Vs + 
{r + 2)^fs + s^/r + e 



< (1 + U)y/~S + Vy/r + | 



□ 



3.2 Sturmian words 

For certain values of r and s, the inequality (J3J holds for all i,f G N even 
when e = 0, and the three coordinates x(t),y(t),z(t) can be determined 
exactly in closed form. In Proposition 13.21 we treat the case r = 2, s = 1. 

By a binary word w = WqWxWi ... we will mean a map N — > {0, 1}; we 
write Wi for the image of i £ N under this map. To each state a £ Rec we 
associate a binary word u> = uv jS (cr), whose n-th term is or 1 accordingly 
as fr S (°~) is /~(.f n ~ :L c) or f + (f n ~ l a). A word u> = w WiW2 . . . is called 
Sturmian if it has the form 

Wn = L(n + l)a + /3J - |na + #] • (10) 
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for real numbers < a, (3 < 1, a irrational. Sturmian words have been 
extensively studied and have many equivalent characterizations; see [3j for a 
survey. The word defined by (fTU|) is called the Sturmian word of slope a and 
intercept [3. 

Proposition 3.2. Suppose r = 2, s = 1, and o = (0,0,0). Then w rtS (a) is 
Sturmian with slope a = \/2 — 1 and intercept (3 = ^ . 

Proof. We have (72,1(0") = and from (J2J), 

x 2 -2y(y- 1) 
Mx,y) = • (11) 

We will show that / n (cr) = o n := (x(n),y(n), z(n)), where 

1 x( n 1 

x(n) = -2|_(n + ~)a|, y(n) = n+— — , z{n) = z (x(n), y{n)). 

Since trivially a n = a and <72,i(0n) = <72,i(c) = 0, by Theorem 12.81 it is 
sufficient to show x{n) < z(n) < y(ri). Note the inequalities 

-2a{n + -) < x(n) < -2a(n+ -) +2; (12) 

tt o; 
(1 _ a ) n _ _ < y (n) < (1 _ a )n - - + 1. (13) 

Since is decreasing on the interval x < 0, we have from (J 12)) 
lV / 1\ . x 2 ,/ l^ 2 



n- ( // + - I -2« U + - + 1< T <^ B + - 
and since 2/(2/ — 1) is increasing on the interval y > 0, we obtain from ()13j) 
(l-a)V- (!-«> + - + — = ((l-a)n--) ((l-a)n---l 



2 4 V 2/ V '2 

< y(y-i) 

< ((l-a)n-|)((l-a)n-| + l) 



1 - a) 2 n 2 + (1 - afn h 



2 
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Also notice that a 2 = 1 — 2a = ~(1 — a) 2 . Now (fTTjl is bounded above by 



2 / 2 1\ (1 — a) 2 2 1 — a 2 a a 2 
z {x,y) < a I n + n + ^J % n + — 2 — " ~ 4 ~~ ~8~ 

= (l-«)n-- + - < y--, 
so z < y as desired. Similarly, (fTTj) is bounded below by 

2 (x, y) > a 2 i^n 2 + n + ^ - 2a + + 1 

(1 — a) 2 9 (1 — a) 2 a a 2 
n l - —n H 



2 2 4 8 

9 7 

= — 2an — a H — > x . 

8 8 

Since x and Zo(x,y) are integers, it follows that z > x. □ 

For general r and s, write w r ^ s = w rjS (0,0,0). It is not true that uv jS is 
Sturmian for every pair r, s. It is a classical theorem of Morse and Hedlund 
(see, e.g., jHj) that a Sturmian word has exactly n + 1 distinct factors 
(subwords) of length n; and it turns out, for example, that w(5, 1) has 70 
factors of length 68. It does not even appear true that every w r , s is eventually 
Sturmian. It does, however, appear that w r<s is Sturmian for a substantial 
number of pairs (r, s). The set of such pairs is quite complex; see figure El 
and the discussion in section 5. 

3.3 Nonlinear bounds and Pythagorean triples 

We now turn to the case r — s — 1 and consider an initial state consisting of 
an interval of occupied sites to the left of the origin. Our next two proposi- 
tions can be seen as rotor-router analogues of the conservation of sx 2 — ry 2 
in the continuum limit of IDLA (see section 2.3). 

Proposition 3.3. Let n be a positive integer and a® = (— n, 0,0). Let o~ t = 

(x(t),y(t),z(t)) = /i,i(cr t -i). Thenx 2 -y 2 < n 2 +n and (x-l) 2 -(y-l) 2 > n 2 
for allt > 0. 
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Figure 2: The lattice path Pi,i(— 4, 0, 0) is bounded between the hyperbolas x 2 ~y 2 = 20 
and (x- l) 2 - (y- l) 2 = 16. 

Proof. For the first inequality, we may assume x 2 > y 2 , i.e. — x > y. Then 
x + 2z < x + 2y < y. Now from lemma 1231 

g-i,i(x, y, z) = x 2 - y 2 - x + y - 2z. 

and since <?i,i(— n, 0, 0) = n 2 + n, at any time t we have 

x — y = n + n + x — y + 2z < n +n. (14) 

For the second inequality, suppose first that x + y < — n. Then 

(x - l) 2 - (y - l) 2 = (x - y)(x + y - 2) > -n(-n - 2) > n 2 . 

On the other hand if x + y > —n, then n + y + 2z > 2,2 — x > x, so from (|14|) , 

(x - l) 2 - (y - l) 2 = x 2 - y 2 - 2x + 2y = n 2 + n - x + y + 2z > n 2 . □ 

Proposition 3.4. Suppose that a,b,n G N are such that a 2 + n 2 = b 2 . In 

the situation of Proposition ^. iA for t = a + b — n, we have x(t) = —b and 
y{t) = a. 
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Proof. By Lemma 12.11 there exists some t such that y(t) = a. Now by the 
first inequality of Proposition 13.31 f° r this value of t we have 

x 2 < n 2 + n + a 2 = b 2 + n < b 2 + b < (b + l) 2 , 

so x > —b. Now by the second inequality of Proposition 13.31 

(x - l) 2 > n 2 + (a - l) 2 = b 2 - 2a + 1 > (b - l) 2 , 

so x < 1 — b. Thus either x = — b or x = 1 — b. In the former case, the proof 
is complete; in the latter case, we will show that x(t + 1) = —b. Indeed, if 
this were not the case, we would have x(t + 1) = 1 — b and y(t + 1) = a + 1, 
but then 

(x(t + 1) - l) 2 - (y(t + 1) - l) 2 = a 2 - b 2 = n 2 , 
contradicting the second inequality of Proposition 13.31 □ 

4 Higher-dimensional analogues 
4.1 The center of mass 

Recall that in dimension two and higher, the rotor-router model depends on 
a cyclic ordering of the set of cardinal directions Ed = {±ei, . . . , ±ed} C 7L d . 
Considering Ed as the set of vertices of a regular octahedron in M. d , if two 
orderings < and <' differ by an octahedral symmetry of E^, the resulting 
rotor-router models will differ by this same symmetry. By applying suitable 
reflections any ordering can be transformed into one satisfying 

(i) ei < e 2 < ■ • ■ < e n ; 

(ii) ei < -ei, z = 1, ... (15) 

For the remainder of this section, we fix an ordering < satisfying (fTB^) (i) and 
(ii). Call a set C C [l,d\ a coclique if the intervals {e e Ed : < e < — e^}, 
% G C are disjoint. 

Theorem 4.1. The center of mass of the set of occupied sites is confined to 
the unit cube {x G M. d : < Xj < 1}. Moreover, if C C [l,d] is a coclique, 
then the center of mass lies below the hyperplane Ylii^c Xi = ^ ■ 
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Proof. For e G E^, let T(e) be the total number of steps taken in the direc- 
tion e by the first n deposited particles. Then - [T(ei) — T(— e,)] is the z-th 
coordinate of the center of mass of the set of occupied sites after n particles 
have been allowed to equilibrate. If the rotor r(x) at the site x G Z d satisfies 
Ci < r(x) < — e», then the site x has ejected one more particle in the direction 
Ci than in the direction —ef, otherwise, x has ejected equally many particles 
in the two directions. Hence 

< T(e<) - T(-ei) <n, 

and dividing by n, we conclude that the center of mass is confined to the 
unit cube. 

If C C [1, d] is a coclique, then every occupied site x satisfies at most one 
of the inequalities < r(x) < — e*, i G C, and so 

< X] t T ( e *) " T (-e0] < n, 

and dividing by n gives the desired inequality. □ 

For example, suppose that d = 3 and < is the ordering e\ < e 2 < —e\ < 
— e 2 < e 3 < — e 3 . Then the sets {1,3} and {2,3} are cocliques, so the center 
of mass is confined to the portion of the unit cube lying below the planes 
x + z = 1 and y + z — 1. 

4.2 Progress toward circularity 

We conjecture that the limiting shape of the rotor-router model, like that of 
IDLA, is a Euclidean ball in M. d . In this section, we prove a much weaker, 
but analogous result, Theorem 14.61 

The discrete Laplacian AF of a function F : Z d — > M is given by 

^F{*) = Y d H F ( x + e )- F ^- 

If AF(x) = 0, then F is said to be harmonic at x. 

Fixing an ordering of write E'd = {ei}^ with ei < ■ • • < e 2 ^. Let 
H m (x) be the total number of times the site x G Z d is been visited by the first 
m deposited particles. The following result shows that H m is approximately 
harmonic away from the origin. 
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Lemma 4.2. I/i ^ 0, then AH m (x) is bounded independent of m and x. 
Specifically, 

~ d+ \~~k- AHm - d+ \ (16) 

Proof. Every time a particle visits the site x, it comes from one of the neigh- 
boring sites x — e, e G E d . When the site x — e, is first visited, the particle 
stays there, and thereafter, the fc-th particle to visit x — e« is routed to x if 
and only if k = i (mod 2d). The total number of routings from x — q to x 
after n particles have been deposited is then at least h [H m (x — e^) — i] and, 
if i < 2d, at most a« = ^ [H m (x — e$) + 2d — 1 — i\. In the case that i = 2d, 
we have a2d = ^ [H m (x — e2d) — 1]? whereas if the site x — e^d has not yet 
been visited, then certainly no routings from x — e<id have taken place, so the 
number of routings is actually a2d+ Summing the contribution from each 
x — €i, we obtain 



2d(2d+ 1) 



Ad 2d 



— ^2 H m {x - e) < H m (x) 



{2d-l){2d-2) 

~\~d ' 2d 



< — — • ^-.Y. "" Ar - () - 



e<EE d 

and this reduces to (fTH|). □ 

There is a unique function G on Z d that is harmonic way from the origin 
and satisfies G(0) = and AG(0) = —1. This G is called the discrete 
harmonic Green's function. Unlike its continuous counterpart, the discrete 
Green's function does not have a simple closed form, and is given in two 
dimensions by an elliptic integral [TS]. However, the discrete Green's function 
does have the same asymptotics as its continous counterpart: in dimension 
3 and higher G is asymptotic to a constant times r 2 ~ d , and in dimension 2 
it is asymptotic to a constant times log^ (cf. [To^ ITS]). 

A first attempt at a proof of the circularity conjecture might run as fol- 
lows. To show that the rescaled set of occupied sites converges to a ball, it 
would be enough to show that the function H m is sufficiently radially sym- 
metric, i.e. to bound \H m (x) — H m (y)\ in terms of ||x|| 2 — ||y|| 2 . Due to 
its asymptotics, Green's function is itself approximately radially symmetric, 
and given Lemma 14.21 we might expect that the function H m should ap- 
proximately coincide with a suitable scaling and translation of G, namely 
H m (0) — AH m (0)G. In two dimensions, however, we immediately encounter 
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the problem that H m is bounded below by zero, while G ~ j- log £ is not 
bounded below. Indeed, the data suggest that G is an excellent approxima- 
tion to H m near r = 0, but a bad approximation when H m is close to zero; 
see figure 0] and the discussion in section 5. Unfortunately, it is precisely 
when H m is close to zero that we need a good approximation. 

We might expect these difficulties to disappear in dimensions three and 
higher, since r 2 ~ d is bounded below by zero. Perhaps surprisingly, however, 
Green's function does not appear to be a very good approximation in the 
three dimensional case, either; see figure |H1 in section 5. 

To avoid the inaccuracies of the Green's function approximations, we will 
take a somewhat different approach. For the sake of simplicity, we treat only 
the two-dimensional model. Fix a map r : Z 2 — > E 2 indicating the rotor direc- 
tion at each point in the plane. Imagine now that several particles are simul- 
taneously deposited at different points in the plane; write v(x) for the number 
of particles at the site x. If a; is a point with v(x) > 1, let x(r, v) denote the 
configuration obtained by routing one particle from x to the neighboring site 
x + r(x), and then changing the direction of the rotor r(x) as dictated by the 
ordering <. A finite sequence of steps (x) = (x±,X2, ■ ■ - Xf.) is said to be ter- 
minating if the configuration (r', z/) = x(r, v) = (x^Xk-i ■ ■ ■ Xi)(r, v) is such 
that v' < 1 everywhere. The following result is a special case of Proposition 
4.1 of [TT], but we include a proof here for the sake of completeness. 

Proposition 4.3. Given a configuration of rotors r : Z 2 — > E 2 , and a map 

Uq : Z 2 — > N indicating the number of particles at each point in the plane, if 
(x) = (x\, . . . , Xk) and (y) = (y\, ■ ■ ■ ,yi) are any two terminating sequences 
of steps, then k = I and the resulting configurations x(tq, Vq) and y(ro, vq) 
are identical. 

Proof. It suffices to show that (x) is a permutation of (y). If this were not the 
case, then reversing the roles of (x) and (y) if necessary, there exists j such 
that the sequence (x 1 ) = (xi, . . . is a permutation of a subsequence of 

(y), but (x") = (xi, . . . , Xj) is not. Then Xj occurs with the same multiplicity 
in (x') and (y), while every x-i ^ Xj occurs with at most the same multiplicity. 
Now setting y(r , v ) = (r, v) and x'(r , v ) = (r', z/), it follows that v\p) < 
z/(p). But (y) is a terminating sequence, so v(p) < 1, and hence p is not a 
legal step after the sequence (xi, . . . , Xj), a contradiction. □ 

Consider now the following procedure. First, m particles are deposited 
simultaneously at the origin, and at each step thereafter, all but one of the 
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particles at each occupied site are routed to neighboring sites, until there 
is at most one particle at each site. By Proposition 14.31 this procedure is 
guaranteed to terminate, and to give the same final configuration (r, v) our 
original model, in which the particles were deposited one by one. Notice that 
at each step in this new procedure, each site ejects approximately equally 
many particles to each of its neighbors; letting H m n (x) be the number of 
particles at the site x after n steps, we see that H mn « if m ri „ 1 + AH m n _i, 
so our procedure has the effect of approximately iterating the operator A+Id. 
Our next two lemmas convert this observation into a precise estimate. 

Lemma 4.4. Let B (x,y) = S x oS y o, where 5 denotes the Kronecker delta. 
Let B n = B n _\ + AB„_i, n > 1. Then 



B n {x,y) =4T 



n 



\(n + x + y) 



n 

\(n + x-y) 



(17) 



where we adopt the convention that 



if k ^ N or k > n. 



Proof. Induct on n. Writing u = \{n + x + y), v = \{n + x — y), we have by 
the inductive hypothesis 

B n _ x {x - l,y) 



1 1— n 



B n _x(x,y- 1) 
B n _i(x,y+ 1) 
B n _i(x + l,y) 



tl-n 



1 1— n 



1 1— n 



n — 1 
u — 1 
n — 1 

u 

n — 1 

u 



n — 1 
v - 1 

n — 1 

v 

n — 1 
v - 1 
n — 1 

v 



whence 
B n (x,y) 



n 
u 



+ 4" 



n — 1 
u 

n 
u 



n 
v 



n 



u 



n — 1 

v 



□ 
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Lemma 4.5. \H m ^ n — mB n \ < 3n. 

Proof. By definition, F m)0 = mS x o5 y o = mB Q . We now show by induction 
on n that 

-3<# TO , n -(A + Id)[# ro)n _i]<3. (18) 

On the n-th step, depending on the direction of the rotor r(z + e) and the 
congruence class of H m ^ n _\{z + e) (mod 4), a point z 6 Z 2 receives between 
| [f/ TOin _i(^ + e) — 4] and | [if m) „_i(z + e) + 2] particles from each neighbor- 
ing site z + e. Moreover, the site 2 ejects all but one of its own particles, 
leaving 

-4 + (A + Id)[ff m , n _i](2) + 1 < F m , n (^) < 2 + (A + Idjfff^n-i]^) + 1, 

and (HH1) follows. □ 

We're now ready to prove the promised weak circularity result. We will 
show that after m particles are deposited and allowed to equilibrate, every 
site in a disc centered at the origin of radius proportional to m 1 / 4 is occupied. 
We will make use of Stirling's inequality 

2n 



and the fact that for any e > there exists N(e) such that for all t > N(e), 



a 



1--J >e' a ~ e (20) 

Theorem 4.6. Let e > 0. In either the type 1 or type 2 rotor-router model 
on Z 2 , if m is taken sufficiently large, then after m particles are deposited 
and allowed to equilibrate, every site in the open disc of radius r = (gl^) 1 ^ 4 
centered at the origin is occupied. 

Proof. Let (x,y) G Z 2 and suppose first that x = y = (mod 2). Write 
\\x + y\, v = \\x — y\. Let r = ^ 'x 2 + y 2 = a/2-u 2 + 2f 2 and let 



u 

2 



a - \r 2 . Note that n is an even integer. By Lemmas 14.41 and 14.51 and 
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Stirling's inequality (|19|). for any m G N we have 



H m Jx,y) > m4 n 



n 



u + n/2 J V v + n/2 



n 



3n 



"" ' r 2 /4 



II II 



4 ^2 

■4 r 



> 



4m 



n 

i=l 



r 2 1 i • u r 2 



-v-l+j 



> 



4m 
4m 



n 



i=i 



4(m+ 1) 
r 2 + M 



3=1 



n 



+ 3 



4(v + l)\_ . 2 
2 



r 2 + 4j 



e 4 r 2 



We now use the inequality (|2Uj) in the cases 
4w(w + 1) 



«2 



4f (u + 1) 



tl = U] 

t 2 = v. 



Let 5 be such that 25 + -rJ^ 
have from (j21j) . 

H m Jx,y) > 



e. For u,v > iV(5) and r < (gf^r) 174 , we 



> 



> 



4m 
e 4 r 2 



exp 



4m 2 + Av 2 + 4u + 4v 



-25 



4m / 2u + 2v\ ., , 
exp 2—^- - -r 



g6+25^2 



4m 



g6+2<5^2 

4m 

g6+25^2 



exp 



exp 



2^2 



3 

2T 2 



sm 



3e^ 



iv(5); 



> 0. 



^ 3 
- — ? 



(22) 



It remains to consider the case that u or v is < N(5). First, let (3 



l-e" 



1/2 



and take m large enough so that mB n (x, y) > 3n for all pairs x, y 
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for which u, v < j3N(5). By Lemma 14.51 H mtn (x,y) > for all such pairs. 
Now suppose one of u, v is > (3N(5). Taking y to — y if necessary, we can 
assume u < N(S), v > /3N(8). In particular, this implies 

Au(u + 1) Au 2 A r 
— < < — = 1 — e 

r 2 v 2 [3 2 

so from (pTjl we obtain 

4m / n Au(u+1)\ ( Av{v + \) ^ 3 2 

4m / 2v 2 2v \ 3 



> -^exp I -4 - 25 - 2 ~" - ) - ~r 2 



Am f 2 \ 3 

> ex P —6 — 2d — —7 rr — -r 
r 2 \ N{5) J 2 



.2 



and we recover equation (J22|) . 

In the cases x = y = 1 (mod 2) and x ^ y (mod 2) the proof is similar, 
taking n = and respectively. □ 

5 Conjectures 

5.1 The Sturmian region 

Figure El shows the set of pairs (r, s) with 1 < r, s < 30 for which the binary 
word w r>s is Sturmian in the first ten million terms. In each of these cases, 
w r>s takes the form (fTU|) with 

We make the obvious conjectures. 
Conjecture 1. 

(i) For (r, s) in the diagonal stripe —4 < r — s < 3, with the single 
exception of the case r = 4, s = 1 the sequence uv jS is Sturmian. 

(ii) If (r, s) is just below the stripe, i.e. r — s = 4, then uv jS is Sturmian 
if and only if r is even. 

(iii) If r — s < —4, then with finitely many exceptions, w r}S is not 
Sturmian. 
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Figure 3: Plot of the set of pairs (r, s), 1 < r, s < 30, for which uv,s) is Sturmian out to 
10 7 places; r increases from left to right and s increases from bottom to top. The square 
of area 1 centered at (r, s) is shaded green (light gray) if w TtS is Sturmian, and blue (dark 
gray) if uv,s is not Sturmian. 



It is likely that (i) and the positive direction of (ii) can be proved in 
the same way as Proposition I3.2| the computations become quite extensive, 
however. Proving the negative direction of (ii), as well as (iii), may be trickier. 
It is not hard to show that if uv, s is Sturmian, then its slope a must be given 
by (j2HJ); indeed, from Theorem 13. II we have — > as t — > oo, and so the 
proportion of u < t for which f(a{u)) = f~(a(u)) approaches 



r 



M 

ry 



1 



M J. y 

r s 



Nil 

ry s 



and we recover ((221) • Likewise, from (jlj) one can deduce that if w ryS is Stur- 
mian, its intercept (3 must be given by (f23|) . 

The fact that w TyS appears to be Sturmian for so many pairs (r, s) suggests 
that an exact description of the triples /™ s (0,0,0) for general r and s may 
be within reach. The subword complexity (number of factors of length n) of 
uv jS would be worth investigating in this connection, since Sturmian words 
are of minimal complexity among aperiodic words (cf. 
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Figure 4: The function H m (r, 0) (middle curve) for the type 1 rotor-router model on Z 2 , 
plotted against the functions F(r) (lower curve) and F(r) (upper curve) of equation Ij24(l . 
Here m = [IOOOOttJ , so that H m has its root near r = 100. 

5.2 Green's function estimates 

The discrete Green's function in two dimensions has the asymptotics (cf. 

EH) 

G(r) = Hogr + | log2 + 7 + 0(-^)^ . 

Set G{r) = ^- (logr + | log 2 + 7). As in section 4.2, let H m (x,y) denote the 
total number of times the point (x, y) G Z 2 is visited by the first m deposited 
particles. As discussed in section 4.2, we might expect H m to coincide closely 
with the function H m (0, 0) — AH m (0, 0)G'. Observe that every time a particle 
visits the origin, it has either been newly placed there, or it has come from 
one of the four adjacent sites (±1, 0), (0, ±1); so by the same argument used 
in the proof of Lemma I4.2| we conclude that AH m (0, 0) is approximately — m. 
Figure |U plots H m (r, 0) against the two functions 

F = H m (0, 0) -mG, F = H m (0, 0) - m \G - . (24) 

Since we constructed F to coincide with H m near zero, and F is unbounded 
below while H m > 0, it is not too surprising that F gives a good approx- 
imation to H m near zero, but a bad approximation near the root of H m . 
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Figure 5: Near its root ro, H m (r, 0) falls off like \r— ro\ x , where A varies only very slightly 
with m. Shown here is H m (upper curve) with m = [f00007rj, alongside the curves H(r) 
and \r- r | A for A = 2.232. 

Interestingly, however, the root of F seems to coincide very closely with that 
of H m . 

Figure El shows H m plotted against H and the function \r — r | A , where 
r ~ 100 is the root of H m and A = 2.232; this latter function gives a very 
good approximation to H m near its root. 

Since Green's function is bounded below on Z n for n > 3, we might expect 
it to give better approximations to H m in higher dimensions than it does in 
dimension two. For the most part, however, these expectations do not seem 
to be bourne out. Figure |B1 shows H m (r, 0, 0) for the three-dimensional rotor- 
router model obtained from the ordering e% < —e% < e2 < — e 2 < e% < —e^, 
plotted against m/r, the principal term of the discrete Green's function in 
three dimensions. 
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Figure 6: The function H m (r, 0,0) for a three-dimensional rotor-router model (darker 
curve), plotted against m/r, the dominant term of Green's function on 1? . Here m = 
|_(4/3)7r • 15 3 J , so that H m has its root near r = 15. 
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